Power.OneWay.Random <- function(n,A,alpha,var.ratio){
 df.1 <- A-1
 df.2 <- A*(n-1)
 F.crit <- qf(1-alpha,df.1,df.2)
 C <- 1 + n * var.ratio
 pow <- 1-pf(F.crit/C,df.1,df.2)
 return(pow)
}